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We perform quantum tomography on one-dimensional polariton condensates, spontaneously oc- 
curring in linear disorder valleys in a CdTe planar microcavity sample. By the use of optical 
interferometric techniques, we determine the first-order coherence function and the amplitude and 
phase of the order parameter of the condensate, providing a full reconstruction of the single particle 
density matrix for the polariton system. The experimental data are used as input to theoretically 
test the consistency of Penrose-Onsager criterion for Bose-Einstein condensation in the framework 
of nonequilibrium polariton condensates. The results confirm the pertinence and validity of the 
criterion for a non equilibrium condensed gas. 
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A gas of bosons, when cooled down to a sufficiently 
low temperature, undergoes a phase transition forming 
a Bose-Einstein condensate (BEC) - a quantum degen- 
erate state in which spontaneous coherence develops be- 
tween the particles PQ. The general definition of BEC 
was proposed by Penrose and Onsager It is based on 
the single particle density matrix p(x,x') — (x)i/}(x')) 
and states that BEC occurs when N c , the largest eigen- 
value of p, is of the order of the total number of parti- 
cles in the quantum fluid. This definition remains mean- 
ingful for inhomogeneous systems and it corresponds to 
the existence of off-diagonal long range order (ODLRO), 
p(x, x') — > n c ^ for \x — x'\ — > oo [3]. In practice 
however, the measurement of the single particle density 
matrix relies on challenging experiments. With ultracold 
atomic gases, the ODLRO was experimentally verified by 
a matter wave interference experiment 4 . To the best of 
our knowledge however, a full experimental reconstruc- 
tion of the density matrix has not been performed yet. In 
this work, we used the favorable properties of microcavity 
polaritons to perform this task. 

Polaritons are quasi-particles that represent the eigen- 
modes of the strong coupling regime between light and 
matter. Such a regime can be achieved in planar semicon- 
ductor microcavities [5], multilayered structures where a 
set of quantum-wells are placed at the antinodes of the 
cavity mode. This configuration allows for efficient cou- 
pling of the photon mode with the exciton resonance, 
yielding polariton cigcnmodcs for the strongly coupled 
system. Polaritons have a bosonic character and favor- 
able properties for Bose-Einstein condensation: an ex- 
tremely light effective mass coming from the photon com- 
ponent and interactions provided by the exciton compo- 
nent. So far, polariton condensation has already been 
demonstrated by several groups [BHS]- 
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the polaritons have a life time, that is in the picosecond 
range. To maintain the polariton density, the losses are 
compensated by a continuous injection of excitons in the 
system. At some threshold injection intensity, the po- 
lariton phase space density becomes of order unity and 
stimulated emission sets in. This leads to a fast increase 
in the polariton density with excitation power and simul- 
taneously, the onset of ODLRO is observed [BJ . While the 
life time of the polaritons is too short for the system to 
reach a global thermal equilibrium, a thermal distribu- 
tion on the lower polariton branch was experimentally 
observed [BJ. 

Thanks to their photonic component, the experimen- 
tal situation for measurements of the single particle den- 
sity matrix is significantly better for the BEC of exciton- 
polaritons. In fact, coherence measurements of polariton 
condensates are much easier than the ones with atomic 
gases. They can be done by optical interference experi- 
ments with the light that is emitted by the microcavity. 
The first unambiguous proof of Bose-Einstein condensa- 
tion of microcavity polaritons was precisely obtained by 
such an experiment [BJ. The long-range spatial coher- 
ence was determined by interfering the light emitted by 
the microcavity with its inverted image, giving access to 
p(x,-x). 

In this Letter we report on the experimental recon- 
struction, through a series of optical interferometric ex- 
periments, of the full single particle density matrix for 
a one-dimensional polariton condensate, i.e. p(x,x') for 
all x and x' . The main motivation for the use of a ID 
polariton state is the significantly smaller size of the den- 
sity matrix with respect to the 2D case. The quantum 
tomography we perform on such a state allows us to test 
the Penrose-Onsager criterion for BEC in the framework 
of Bose-degenerate gas of polariton quasi-particles. We 
provide convincing evidence of the pertinence and valid- 
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ity of the Penrose- Onsager criterion even in the case of 
non-equilibrium BEC, as it is the case for the polariton 
system. 

The sample is the same CdTe planar semiconductor 
microcavity of our previous works [5] . The sample fea- 
tures a pronounced disorder potential, naturally arising 
as a result of the epitaxial growth process: particular 
linear disorder landscapes have already been proved to 
allow for single-energy ID polariton condensates [TP] , 
The quantum tomography - the density matrix recon- 
struction - requires the determination of both the am- 
plitude and the phase of the first-order coherence func- 
tion g^(x,d) = p(x,x + d) of the condensate, where x 
represents the spatial coordinate along the ID conden- 
sate direction and d is a given spatial separation between 
points of the condensate. Along with the coherence mea- 
surements, the condensate order parameter itself requires 
a characterization in terms of its density and phase struc- 
ture. In order to measure the required quantities we make 
use of advanced interferometric techniques. 

First we proceed to the determination of the first order 
coherence function. A schematic depiction of the exper- 
imental setup used for this purpose is shown in Fig. [T] 
The setup consists of a photoluminesccncc experimen- 
tal apparatus in reflection configuration. The sample is 
held in a He-cooled optical cryostat at approximately 4 
K. We shine a non-resonant Ti:Sapphire quasi-CW laser 
through a 0.5 NA microscope objective on the surface of 
the sample to create a polariton population in the sys- 
tem. The non-resonant excitation scheme ensures that 
no coherence is imprinted in the system by the laser onto 
the polariton gas. The PL signal emitted by the sam- 
ple is collected through the same microscope objective 
and is send to an actively-stabilized Michclson interfer- 
ometer [5J, which represents the core of the detection 
scheme and allows for high-stability phase measurements 
required for the purpose of this work. The interferom- 
eter, differently from our previous work on ID conden- 
sates |l(Jj . is in the mirror-mirror configuration that is 
most convenient for the purposes of quantum tomogra- 
phy, as detailed in the following. The real-space conden- 
sate density is imaged by using a CCD camera placed at 
the output of the interferometer. 

Starting from low excitation power, we observe emis- 
sion from a classical gas of polaritons. By increasing 
the pump power, once we cross the threshold power of 
35 \iWjcm} , we observe the spontaneous onset of conden- 
sation, characterized by the emission of a macroscopically 
populated single-energy polariton state. By slightly tilt- 
ing one mirror of the interferometer (displacement con- 
trol mirror in Fig. [TJ we are able to set the overlap condi- 
tions between the two arms with micrometer resolution. 
Once an overlap condition is set, corresponding to a spe- 
cific value of the d parameter in the density matrix, a scan 
over 6ir of the delay between the two arms of the interfer- 
ometer is performed by shifting the other mirror of the 
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FIG. 1: (a) Experimental setup. The non-resonant excita- 
tion laser is focused through a 0.5 NA microscope objective 
on the sample (hold in the cryostat at ~ 4K). The same ob- 
jectives collects the PL signal emitted by the sample, which is 
then sent at the input to the Michelson interferometer in the 
mirror-mirror configuration. One mirror is actively stabilized 
and controls the delay between the two arms; the other mirror 
is adjusted to realize the desired displacement d between the 
two arms. Polariton density below (b) and above (c) conden- 
sation threshold. Condensation is observed to occur in a ID 
state. 

interferometer (delay control mirror in Fig.[l]). The scan 
of the delay is computer controlled by taking advantage 
of the active stabilization mechanism of the interferome- 
ter [5] . This experimental procedure yields a discrete set 
of interferometric images, one for each delay condition of 
the interferometer arms. Thus, for each spatial point of 
the condensate, one obtains an intensity modulation as 
a function of delay. By fitting such intensity modulation 
we retrieve a map of the amplitude of the g^(x,d). 

We consider a finite set of values for the displacement 
d, corresponding to steps of « 1.2 /im, covering the whole 
length of the ID condensate. The measurements are re- 
peated for different excitation powers, below and above 
threshold. In particular, above threshold, we made sure 
to have only a single-energy condensed state emitting and 
to be far from the onset of multimode condensation. 

As previously mentioned, access to the phase infor- 
mation is also required in order to fully reconstruct the 
density matrix. In order to experimentally assess this 
quantity, under the same excitation conditions, we make 
use of a different detection scheme. Such scheme basi- 
cally consists of a Mach-Zehnder interferometer in which 
one arm is the polariton condensate density and the 
other arm is a magnified version of it. A small part 
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of the ID condensate, enlarged by four times, acts as 
an approximately flat phase reference and is overlapped 
with the whole condensate density, similarly to what 
was done in Ref. From the resulting interferogram, 
the spatial variation of the phase of the polariton field 
Acf>(x) — 4>{x) — 4>(xq) is extracted. The density matrix 
is then constructed from the absolute value and phase 
variation as p{x,x') = ^{x.x')^^'^^^ . The fact 
that a non-constant phase exists for the single parti- 
cle density matrix is due to the nonequilibrium state 
of the polariton quantum fluid, that breaks time invari- 
ance. Consequently, in contrast to equilibrium conden- 
sates [12], the steady state does not correspond to the 
ground state of the system and its wave function is not 
guaranteed to be real. Within the polariton condensate, 
particles flow proportional to the gradient of the phase 
are allowed. This peculiarity of the nonequilibrium po- 
lariton condensates has made it possible to observe quan- 
tized vortices without setting the system into rotation 

nam]. 

The density matrices obtained through quantum to- 
mography are shown in Fig. [2] Each column of the ma- 
trices represents the value of coherence along the con- 
densate, between spatial points that are separated by a 
distance d! ~ d- 1.2 fim (d = corresponds to full overlap 
whilst for d = 18 the two arms have a negligible over- 
lap region). The amplitude of the g^(x,d), for a pump 
power below (20 pW/cm 2 ) and above (60 ijW/cm 2 ) con- 
densation threshold, is shown in Fig. |2ja) and (b) re- 
spectively. The difference between the two situations 
is striking: in the classical gas the coherence drops to 
zero within a few microns of displacement and becomes 
completely negligible already at step d — 4; in the con- 
densate the off-diagonal long-range order (ODLRO) is 
well established and the coherence has significant values 
over the whole length of the condensate, even up to step 
d = 17. In Fig. |2jc) and (d), below and above condensa- 
tion threshold, the phase of the density matrix is shown. 
Below threshold the relative phase is not well defined and 
it is seen to randomly fluctuate already for very small 
displacement values. On the contrary, above threshold 
a very well defined phase structure is identified, with a 
phase profile that varies smoothly as a function of spatial 
position x and displacement d. 

The correct phase of the density matrix is important to 
obtain the correspondence between the Fourier transform 
of the density matrix and the momentum distribution of 
the polariton gas. 
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n(k) = / dxdx' e lk{x ~ x '^ p(x,x'). 



(1) 



The momentum distribution reconstructed from the den- 
sity matrix is shown in Fig. [3] When the phase informa- 
tion of the single particle density matrix is retained, the 
momentum distribution is correctly peaked at a negative 
momentum. Instead, when the phase information is not 
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FIG. 2: Density matrix reconstruction for the polariton sys- 
tem: amplitude of the g^(x,d) in the case below (a) and 
above (b) condensation threshold, respectively for pump pow- 
ers of 20 pW/cm 2 and 60 fiW/cm 2 . The corresponding phase 
structure is shown in (c) and (d) . A high degree of long-range 
order is present in the condensed phase as well as a defined 
phase structure. Below threshold the measured data stop at 
step d = 7. 



included and a constant phase is assumed all over the 
condensate, the calculated momentum distribution ap- 
pears centered around zero. The shift of the peak of the 
momentum distribution toward negative momenta corre- 
sponds to the experimentally measured momentum dis- 
tribution and it is consistent with the measured phase 
gradient along the one-dimensional condensate that cor- 
responds to a flow of polaritons with a finite wavevec- 
tor. Comparison with the directly measured momentum 
distribution (see Fig. [3jb), dashed thick line) however 
shows that the reconstruction of the momentum distri- 
bution from the density matrix is far from perfect: the 
experimental momentum distribution shows much more 
structure than the reconstructed one and its overall shape 
is much broader. We suspect that this is due to the fi- 
nite spatial resolution in the phase measurement. This 
results in an underestimation of the spatial phase varia- 
tions and consequently a momentum distribution that is 
too narrow. 

The diagonalization of the measured density matrix 
yields a set of eigenvalues, representing the population of 
the corresponding eigenstates for the system. The cal- 
culation of the eigenvalues allows us to directly validate 
the Penrose-Onsager criterion for BEC. For a classical 
gas, the quasi-particles are distributed among the differ- 
ent states available in the system, corresponding to a set 
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FIG. 3: The momentum distribution computed from the 
single particle density matrix containing the phase informa- 
tion (full line) and from \rho\ (dashed line), below (a) and 
above (b) threshold. When the phase information is omitted 
[Th. (np)], the momentum distribution is incorrectly cen- 
tered around zero. With phase information [Th. (wp)], the 
distribution peak is shifted towards the left, towards the same 
direction as the measured k-space distribution (thick dashed 
line) . 



of comparable non-zero eigenvalues. This is indeed what 
we find from the experimental data. Fig. |3Ja) shows the 
occupations of the single particle states. As expected 
it has the shape of the Maxwell-Boltzmann distribution 
below threshold, corresponding to a Gaussian correlation 
function in real-space, as shown in Fig. |4|c). The densi- 
ties of the three most occupied states are also plotted in 
Fig. |4jc). They show the spatial dependence of standing 
waves with increasing number of nodes. 

On the contrary, upon polariton condensation in a 
single-energy state, only one state is macroscopically 
populated so that one eigenvalue should be found signifi- 
cantly different from zero and dominant over all the oth- 
ers [Fig. Qb)]. This implies that a single-energy conden- 
sate was created. Co-existence of multiple condensates 
on the other hand would lead to several large eigenvalues 
[T5], [TH]. The position on the sample was actually care- 
fully chosen for a single condensed state to occur. When 
scanning the position of the excitation spot, a place was 
chosen where a single line was dominant in the frequency 
spectrum of the condensate. The present analysis of the 
eigenvalues of the density matrix shows that this cor- 
responds to a single-energy condensate. Indeed, a very 
good agreement between the macroscopically populated 
eigenstate and the measured polariton emission is found, 
as shown in in Fig. |4|d). 

Unfortunately, we find non-physical negative occupa- 
tion for the last few eigenvalues of the density matrix 
spectrum [see Fig. |4jd)]. We understand this spurious 
effect as a result of the experimental error on the deter- 
mination of the coherence: such error is too large to at- 
tribute a physical meaning to the states that correspond 
to the weakly occupied states. The contribution of the 
condensate to the density matrix p c = N c <p* (xi)<j){x2) 
is so large that the experimental error on the remainder 
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FIG. 4: Results from the diagonalization of the density ma- 
trix. Below threshold (a,c): eigenvalues (a) and the eigen- 
vectors (c) of the three states with the highest occupation. 
Above threshold (b,d): the eigenvalues (b) and the state of 
the condensate (d). 



p — p c is too large to yield reliable results. 

In conclusion, in this work we have performed a full 
measurement of the single particle density matrix of a 
nonequilibrium polariton condensate and we have com- 
pletely determined the complex order parameter. The 
consistency of the Penrose-Onsager criterion for BEC 
above the threshold laser power was experimentally ver- 
ified. This results give rigorous support and strengthens 
the analogy with the cold atoms BEC, reaffirming once 
more the central role of polaritons in the study of con- 
densates and related phenomenology. 
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